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Error probability is considered for binary signaling through a midtipath 
channel in which (i) the receiver observes a waveform comprising white 
Gaussian noise and the sum of (perhaps several) time-delayed, frequency- 
shifted, Rayleigh-faded versions of the transmitted waveform, (ii) the 
receiver decides with minimum error probability which of the two possible 
transmissions was sent. Results given herein for the exact minimum error 
probability necessarily depend upon a number of parameters and are 
cumbersome to use. By introducing bounds on the error probability, de- 
pending upon bounds on spectra of certain matrices, the number of param- 
eters is reduced and the less cumbersome residts become applicable to any 
one of a set of channels rather than to just one channel. The error-prob- 
ability bounds are presented in terms of values of the distribution function, 
derived herein, of the difference of two chi-square random variables. The 
bounds are sharp when the spectra are narrow. For the case of widely 
orthogonal signals, any version of one possible transmission being orthogonal 
to any version of the other tra/ismission, the bounds are given as a set 
of universal curves plotted versus signal-to-noise ratio for various values 
of the number of paths and of the spectral width of certain matrices. Spectral 
bounds can easily be computed when the versions for each transmission 
are nearly orthogonal. Returning to the general case, another bound is 
derived, by a technique due to Chernoff, which does not explicitly require 
spectral bounds which may neither be readily available nor be accurate 
approximations of eigenvalues. This bound is not as sharp as the previous 
bound for the case of small spectral width, but has promise for the large- 
width case. 

I. INTRODUCTION 

This paper considers error probability for the optimum reception 
of binary signals transmitted through a multipath channel having 
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P paths.* One of two possible signals is transmitted; the received 
waveform is the sum of P Rayleigh-faded, time-delayed, frequency- 
shifted versions of the transmitted s igna l, plus white Gaussian noise. 
That is to say, if the complex signal \^2E m x m (t) is transmitted, m = 1,2, 
the contribution to the received waveform from the pth path is 

yM)p) = V2K, a p x„,(t - t v ) exp [i(2*f v t + *>,)], 
where a p , cp p , r„ and f p are the Rayleigh-distributed amplitude, the 
uniformly-distributed phase, the fixed time delay, and the fixed fre- 
quency shift associated with the pth path. The received waveform is 

zM) = JLyM;v)+n{t), 

p=i 
where n(t) is white Gaussian noise. 

The above multipath situation is a special case of a more general 
communications situation in which a receiver observes a sample z{t) 
of a zero-mean complex Gaussian process on the time interval [0, T], 
the covariance function (z{s)z*(t)),„ having been selected from a set 
of two distinct functions by chance according to the prior probabil- 
ities {«„,}, m = 1, 2, and the other second-moment function (z(s)z(t)) m 
being zero. The receiver is to be designed so that its decision 
upon one of the two possible hypotheses is made with mini- 
mum average error probability P e , where P c = 2a,„P c (ra) and 
P e (m) is the probability, when covariance indexed in is true, of 
deciding otherwise. 

The receiver-design problem has been treated in Ref. 1, rigorously 
demonstrating that optimum processing involves quadratic filtering. 
However, the filter kernels, being the solutions of integral equations, 
are difficult to determine in general; moreover, the error probability 
is not evaluated. For the multipath channel, the first difficulty is over- 
come in Ref. 2 and the evaluation of binary error probability is 
considered in the present paper. 

Section II presents the theory of a method that can be used to 
calculate error probability exactly. However, it is quickly appreciated 
that error probability depends in a cumbersome fashion upon a 
large number of parameters including the path strengths and the 
scalar products of the versions. To simplify this situation, this paper 
introduces bounds on the error probability which depend upon bounds 
on the spectra of certain matrices, the eigenvalues of which determine 



* Each path could comprise a multitude of randomly phased subpaths having 
essentially the same delay and frequency-shift parameters. 
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error probability exactly. Thus, the bounds are applicable to any one 
of a set of channels rather than to just one channel. 

Section III presents these error-probability bounds in terms of 
values of the distribution function of the difference of two chi-squarc 
random variables and then derives this distribution function. More 
specific results are obtained in Section IV for the case of widely 
orthogonal signals, any path's version of one of the two possible trans- 
mitted waveforms being orthogonal to any path's version of the other 
waveform. Here, easily computed spectral bounds can be given for 
the case in which the versions under each hypothesis are nearly 
orthogonal. Section V considers the case of well-resolved paths, making 
contact with diversity theory (Ref. 3, Chap. 7), and the case of on-off 
keying. 

The error-probability bounds considered above require spectral 
bounds which may not always be easily computed and which may not 
be accurate approximations of the eigenvalues. A bound that circum- 
vents these difficulties is obtained in Section VI with a technique due 
to Chernoff. Comparison of this bound with previous bounds is car- 
ried out analytically only for the case of well-resolved paths, but 
qualitative comparison is made for more general cases. 

IT. PROCEDURE TO OBTAIN ERROR PROBABILITY IN THE GENERAL CARE 

2.1 Notation 

The binary situation is a specialization of the case of M-ary 
signaling through the multipath channel in which the received process 
z(t) can have one of M possible covariance functions, (z(s)z*(t)) lll} 
m = 1, 2, . . . M, of the form 

r 

2E m £ aMs, m)b*(t, m) + N 8(s - t), 

a degenerate kernel plus a white-noise kernel (Ref. 2). Here b p (t, m) = 
exp (i2irj p t)x,„(t — t„) is a time-doppler-shifted normalized version of 
the transmitted signal \/2E m x m (t); the path with index p has an 
average cross section of a v units, a delay of r p seconds, and a doppler- 
shift of /„ Hz. We put 

/ dt |.r„,(/) | 2 = f dt I Kit, m) | 2 = 1, 

so that the average energy received from the medium is 
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| f dt 2E m 2 a v | b p (t, m) | 8 - E m Z *v - *- , 

since we put 23 o-„ = 1. 

The above covariance function can be written 

2# m 5(s, m)<rb*(t, m) + JV 5(s - t), 

where b(t, m) is a vector with pth component b v (t, m) and a is the diagonal 
matrix with pth entry <r v , with tr a = 1. 

The optimum receiver decides according to the value of m that 
corresponds to the largest of M test statistics computed as follows. 
For each value of ra, the receiver first generates the column vector 
Z(m) = No* / dt z(t)b*(t, m) and then evaluates a test statistic com- 
prising a Hermitian form in Z(m) plus a bias constant. This test statistic 
is 

[{N Q /2EtfZ{m)]X2EJN Q )H(m)[(N /2E$Z{m)-\ + {N /2E x )6(m), 

where the Hermitian combining matrix is 

(2E m /N )H(m) = (2EJN )[(2E m /N )B(m) + (r -1 ]' 1 , 

the bias is given by, 

.,. _ a m det [(2E 1 /N )B(1) + a' 1 ] 
6{m) - log^ det [{2E j NQ)B{m) + ,-»] . 

5(m) is the correlation-function matrix / ^ &*(i, m)h(t, m), and the 
hypotheses are ordered so that E x = max E m . The above test statistic 
is obtained from that given in Ref. 2 by subtracting log [a t det -1 a 
det -1 # -1 (l)] and multiplying all resulting terms by N /2E 1 . 

The above test statistic has a certain intuitive appeal. The components 
of the vector Z(m) are the correlations of the received signal against 
the noise-free versions of the transmitted signal that would occur 
when message m is sent. That is to say, Z(m) provides a measure 
of the projection of z(t) on the P-dimensional subspace spanned by 
these versions. Moreover, the test statistic is a measure of the likelihood 
that this P-dimensional subspace is in fact the correct subspace. Then 
the optimum receiver strategy is decision according to the most likely 
of the M possible subspaces. Also, since P dimensions are involved, 
it might be anticipated that the results are related to the case of P-fold 
diversity, cf. Section 5.1. 

Henceforth, only the binary case, M = 2, is considered. In this 
case, decision according to the larger of two test statistics is equivalent 
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to decision according to the sign of their difference. The decision events 
can then be written in terms of one Hermitian form in a composite 
Gaussian vector 



as follows. Let 



Q = 



2E t 



//(I) 



N 

OpxP 



Opxp 



where O^x^ is the p X p zero matrix. Then the receiver decides upon 
m = 2 when Z f QZ is less than (N /2E 1 )6(2), and decides upon m = 1 
otherwise. 

The conditional error probabilities are thus 

P.(l) = Pr [Z'QZ < (N /2E)6 | 1} = F t (f|)*] , 

P.(2) = Pr {Z'QZ > (N /2E)6 | 2} = 1 - F 2 (^jdj , 

where E = E x , 8 = 6(2), and F m (x) is the distribution function of 
Z*QZ conditioned upon the with hypothesis. 

2.2 The Fundamental Matrices 

Since Z^QZ is a function of a Gaussian vector, the distribution 
function F„(x) is determined by the conditional mean, (Z) m , which is 
the zero vector, and by the conditional co variance L(m) = (ZZ i ) m , 
the other second-moment matrix (ZZ) m being the 2P X 2P zero matrix. 

The conditional covariance matrix L(tn) is evaluated as follows. Let 



L(m) = 



L l \m) 
[L 2 \m) 



V\m) 
L 2 \m)\ 



where L ik (m) = (NJ2E l )(Z(j)7J(k)) m . Then, by the definition of Z(j) 
and interchange of operations, we obtain 



(Z(j)z\k)) m = yJ! ds dt b *(*> JMs)z*(t))M, k) 
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= W II ds dt 6 * (s ' fi 2Em ^ s ' m ) ff6 *(*> m )% k ) 

+ jf ds dt b*(s, j) 8(s - t)h(t, ft) 

OBJ 

= J^B(j, m)aB(m, k) + B(j, k), 

where B(j, in) = / ds b*(s, j)b(s, in) is a cross-correlation matrix. 
Hence, 

L*\m) = (EJEOBd, m)*B(m, k) + (No^EJBU, k). 

Similarly, it is found that (ZZ) m is the 2P X 2P zero matrix. 
For future computations, it is convenient to write 

^ \ Q 22 J ' 
where 

Q 11 = Br\l){I + (No^EJlBil)*]- 1 }- 1 , 

Q 22 = -(E 2 /E 1 )B-\2){{E 2 /E 1 )I + (N /2E l )[B(2)a]- 1 \- 1 . 

2.3 The Characteristic-Function Method 

To obtain the distribution, consider the conditional characteristic 
function 

<p m (t) = <exp (itZ'QZ)) m . 

It is well known, e.g., Ref. 4, that 

*>„(<) = det -1 [/ - itL(m)Q] = Yl [1 - ^X,(m)] _1 , 

k 

where {A*(m)| is the set of eigenvalues of the matrix L(m)Q. The 
eigenvalues are real, since L(m)Q is similar to the Hermitian matrix 
l)(m)QL h (m). 

The distribution function can now be obtained from the characteristic 
function. As a preliminary, it is noted that the characteristic function 
(1 — it\)~ n corresponds to one of two distribution functions, according 
to the sign of X. When X is positive, the distribution function is 



/: 



dx m «- «P (-»A) . J/toA, n - 1) to > 0) , 
x(n - 1)! 1 to<0), 

= U(y)I(y/\,n - 1), 
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where U(x) is the unit step function (unity for x > 0, zero for x < 0, 
one-half for x = 0) and where 

J (y> n ) = —, / dx x " e ~ x = ! - e ~' S IT, 

is the incomplete gamma function. Similarly, when X is negative, the 
distribution function is 



/: 



dxU(-x) ( - X ' y ~ ,eXp(x]Xr) 
ax uk x) |x , n(n _ 1}! 




(V > 0), 

(y < o), 



= 1 - U(-y)T(ij/\,n - 1). 

To obtain the distribution function of Z^QZ, the characteristic func- 
tion is expanded into its partial fractions. Each term will be propor- 
tional to (1 — it\y n for some n, and corresponds to a term in the expan- 
sion of the distribution function. For example, when all eigenvalues 
are distinct, the expansion of the characteristic function is 

d k (m) 



where 



* W ~ f 1 - tfl^m) ' 



The expansion of the distribution function F m (x) is then 
£ d k (m)U(x)I(y/\ k (w),0) 

It :Xt(m) >0I 

+ E d k (m)[l - U(-x)I(ij/Um),0] 

[k :\t(m) <0| 

In the case of a degenerate spectrum, an eigenvalue X with multiplicity 
r contributes the sum ^Z«=i -Ml ~ ^X)~" to the expansion of the 
characteristic function, and the corresponding part of the distribution 
function involves /(-, n) for n = 0, 1, 2, • • • , r — 1. 

It should be observed that the general approach of summing distribu- 
tion functions corresponding to partial fractions is fully equivalent to 
inverting the characteristic function by contour integration, the ap- 
proach used by Turin 5 for a similar problem. (When all poles are simple, 
the expansion coefficients {d k (m)\ are residues of the poles.) 



1608 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1967 

III. UPPER AND LOWER BOUNDS ON THE ERROR PROBABILITY 

3.1 Error-probability Bounds from Degenerate-spectrum Variables 

Exact computation of error probability involves considerable nu- 
merical work in computing eigenvalues followed by evaluation of cum- 
bersome formulas. Moreover, an often inordinately large number of 
independent parameters must be specified. To simplify this situation, 
we consider bounds on the spectrum of L(m)Q rather than the spectrum 
itself. With a technique suggested in Ref. 6, we can obtain error- 
probability bounds. Although we do not obtain the error probability 
itself, the error-probability bounds apply to not just one channel but 
rather to any channel for which the spectral bounds are met. 

Observe that the characteristic function is precisely specified by 
the spectrum of L(m)Q. This spectrum is the same as the spectrum 
of / diag [Xi(m), ■ ■ , \ 2 p(m)], where I plays the role of a covariance 
matrix and the diagonal matrix plays the role of a matrix of a Hermitian 
form. Hence, the distribution of Z^QZ is the same as the distribution of 

ip 
q(m) = 2 **(»») | 2* T, 

where \z k ) are complex zero-mean Gaussian variates with covariance 
matrix (z s z%) = S ik , (z,z k ) being zero. 

Suppose bounds on the eigenvalues are available. That is to say, 
suppose it is known that the positive eigenvalues satisfy 

M ^ X*(to) ^ H, (la) 

and that the negative eigenvalues satisfy 

-v ^ \ k (m) g -v, (lb) 

where the ^'s and v's are positive numbers that depend on m. Then, 
a lower bound on q(m) is the degenerate-spectrum random variable 
q(m), defined by 

P 2P 

q{m) = n J2 I s* I 2 — 5 £ I ** I 2 - 

t=l k=P+l 

Note that we have used the fact that the number of positive eigenvalues 
and the number of negative eigenvalues are the same, see Appendix A. 
Similarly, an upper bound on q(?n) is provided by the random variable 



q(m) = p. 2 I 2* |* "~ " S 



2P 

' z. 
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Since q(m) ^ q(m) ^ q(m), it follows that 

Pr ($(m) ^ y) < F,„(y) = Pr {g(m) fg y) ^ Pr \q{m) < y). 

Evaluation of these bounds requires the distribution function 
G(y; P, a) of the degenerate-spectrum random variable 



£ | Z k | 2 — a 2 



2* 



which is the difference of two clii-square variables each with an even 
number of degrees of freedom. The bounds become 

where we use ?/ = (N /2E)d and reiterate that the m's and j/'s depend 
on m. 

It is anticipated that these bounds are sharp when the spectrum 
is narrow, the spread of the positive spectrum being much less than 
any positive eigenvalue and similarly for the negative spectrum. Also, 
when itself is not precisely known, but bounds ^ ^ are available, 
the distribution function is bounded by 

G[{nY'y; P, .(m)" 1 ] ^ F m (y) g (?[(m)"V; P, Km) -1 ], (2) 

where?/ = (NJ2E)6 and £ = (N /2E)6. 

3.2 Distribution of a Degenerate-Spectrum Variable 
It will be demonstrated that G(y; P, a) equals 
P - 1 + k\( l 



1 +aJ S V * " ! " « 

when y < 0, and equals 

fp - 1 + feV i 



1 - Jl-L«-I,p- 1-fcl (3a) 



when y > 0. 

Before doing so, note that when y < 0, the parameter a serves as 
a scale size for y in the argument of I(x, n), but that this is not true 
when y > 0. Nevertheless, a does act as a scale size in the following 
way. A power-series expansion of I(x, n) yields 



k 



I(y, P - 1 - k) 



(1 + «)* 

1 + J \J (P - k)\ f=L n\ P-k+n J 
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and when a«l, the factor (y/a) p ~ k determines the small-?/ behavior. 
Also, this result exhibits [a/ (I -f- a)] p as a factor for the case y > 0, 
in agreement with the expression for y < 0. 

To find G(y; P, a), we consider its characteristic function 

(1 - ^)- p (l + ita)- p . 

Let the partial-fraction expansion of this characteristic function be 

m-0 n=0 

To evaluate A P _ m , multiply by (1 — it) p and let 1 — it = t to obtain 
(1 + a - ar)- p = £ A P . m r m + r P £ B P . n {l + a - ar)-'^ 1 . 

m = n-0 

Since the second sum is analytic at t = 0, we have exhibited the Taylor 
expansion with remainder. But 



(1+a-ar)-' - (1+^1 ~ fT^ T 



-j» 



1 V f; (P + * - lVr*-)V 



J + a/ i= \ fc / M + a 

where we have used (7) on page 2 of Ref. 7. Hence, 

Ap - = [i~+ij [ m XnrJ • 

Similarly, to obtain B P _„ , multiply by (1 + ita) and let 1 + ita = r 
to obtain 

a al ~ " Oti 

Reasoning as before, it is seen that 

a YfP -\ - I 



B P -„ — i . * \ i \ 

J + al \ n Al + al 

Collecting these results, it is seen that the characteristic function is 



' _I (P + k- l\( 1 Xir/ 



E 



a \ n , .. >.-(/>-*) 



(1 + tfa)' 



i_u \ ^ / \1 + a/ _\1 + a 

a 



This immediately establishes the distribution function (?(?/; P, a). 
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IV. WIDELY ORTHOGONAL SIGNALS 

4.1 Matrices for the Two Hypotheses 

We consider the special case in which the signals are widely orthogonal, 
B(l, 2) = B{2, 1) = Op X p . That is to say, all time-doppler shifted 
versions of one signal are orthogonal to all such versions of the other 
signal, a situation that would prevail in frequency-shift keying with 
widely separated frequencies. In this case, 

V\m) = 8,*[frBG\ m)aB(m, k) T ^B0', *)]• 

The "diagonal" form of the covariance matrix L{m) and of the matrix Q 
implies that the spectrum of L(m)Q comprises the specturm of L ll (m)Q u 
together with the spectrum of L"(m)Q". This can be seen by employing 
the formulas of Schur (Ref. S, pp. 45-46) to reduce the determinantal 
equation det [L(m)Q — XT'] = from order 2P to order P. For m = 1, 

L u (l)Q u = B(l)<r, 

L 22 (1)Q 22 = -(No/2E i )(E 3 /E l ){(E i /E l )I + (N /2E 1 )[B(2)aT i r . 

For m = 2, 

L U (2)Q U = (N /2E t ){I + (No^EMBa)*}- 1 }- 1 , 

ir-(2)(? 2 = -(E 2 /E 1 )B(2)<r. 

It should be observed that the spectra of the above matrices are 
simply related to the spectra of B(\)a and of B(2)a. When E 2 = E^ = E, 
the spectrum of L"(l)Q" is { - (N /2E) (1 + (JVo^W; 1 )" 1 ), where 
{8 k } is the spectrum of B(2)a. Similarly, the spectrum of L 11 (2)Q 11 
is {(N a /2E)(l + (No^E)^; 1 )- 1 ], where {u k \ is the spectrum of B(\)a. 

Second, it should be observed that when E 2 = E t = E, the forms 
of the matrices for the cases m = 1 and m = 2 are the same, with the 
roles of positive and negative matrices interchanged. To compute 
error probability for m = 1, we use the distribution function Fj(x); 
for m = 2, we use the conjugate distribution 1 — F 2 (x) which can be 
expressed as P\—Z^QZ < — x | 2}, the distribution function of the 
negative of the original variable evaluated at — x. Introduction of this 
random variable for the case m = 2 reverses the roles of positive and 
negative matrices, the net effect being that for both m = 1 and m = 2 
the positive and negative matrices have the same forms. 
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4 .2 Bounds on Spectra, 6, and Error Probability 

It is clear that spectral bounds on L(1)Q can be obtained from spectral 
bounds on B(m)<r, m = 1, 2, and similarly for L(2)Q. Consider the 
bounds on L(1)Q when E x = E 2 = E. The positive spectrum is bounded 
as follows: 

n = co ^ min u k ^ X(l) ^ max u k ^ u = /x, 

and the negative spectrum is bounded as follows : 

-v = -{N /2E)[l + (N /2E)(^- l y 1 ^ X(l) 

X(l) g -(N /2E)[1 + (N /2E)(S)~ l T i " -fi 

where 3 = rnin 5* ^ max 5* ^ 5. 

Moreover, bounds on can also be obtained. When E x = E 2 = E 
and a x = a 2 = § (equilikely signals), 

Since a determinant is the product of the eigenvalues of the matrix, 
we have 

(M * = (Mio g n^|f§- 

Thus, an upper bound is 

(N /2E)6 = (N (I /2E)P log ^ffij Uy , 
and a lower bound is 

( No nm = (iv./2E)p ^i \X%1m ' 

Recall that the distribution function F 1 [(N /2E)6] is bounded from 
above by G[([i)~\No/2E)d; P, vfe)' 1 ]. Further, suppose that the 
spectra of B(l)<r and B(2)a are narrow about the nominal value 
(1/P) tr B(m)ff = (1/P) tr <r = (1/P). We can put 

- 5 1 ±J * _ 1 ~ g m 

cu = 5 = — p — , u - 5 — p , W 

where /3 is the fractional spectral half width. Then, the parameters 
required to compute the upper bound on the distribution function are 
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(»)-\N /2E)6 = j-L- (N P/2E)P log Itg + gggg ■ (-) 



Ka)' 1 = rzi ( Ar oP/2£) 



1 + r+^ (Ar ° p/2£) ] *■ (5b) 



Similarly, the distribution function F 1 [(N /2E)8] is bounded from below 
by G[{nY l (N /2E)6; P, wfr)~% The parameters required for this bound 



are 



(,r(N /2E)6 = ^ (WV^P bg ^g + gg gj , (5c) 



zOO"' = rT/3 (iVoP/2£) 



1 + r=i w/ 2 ^ 



• (5d) 



Having considered the case m = 1, the bounds for the case m = 2 
are apparent. Considering the random variable —Z 1 QZ with assumed 
known, the positive and negative spectral bounds are precisely the 
same as for the case m = 1, and the upper bound is 

G[-{^)-\NJ2E)e;P,v{tt)~ x ] 

whereas the lower bound is G[-{p.)~\N o /2E)0; P, rfr)' 1 ]. But 6 is 
unknown, and the upper bound is given by replacing — 6 by 6, and the 
same result is obtained as previously; similarly, the lower bound is 
given by replacing — 6 by 6. In short, the bounds apply to both cases, 
m = 1 and 2. 

The numerical values of these bounds are given in Figs. 1 to 3 as 
functions of 2E/N P (the signal-to-noise ratio per path) for various 
fixed values of £ (the fractional spectral half-width) and P (the number 
of paths). The curves are nested with respect to values of the fractional 
spectral half-width /3; an increase of /3 always yields an increase of the 
upper bound and a decrease of the lower bound. A measure of the sharp- 
ness of the bounds (given a nominal value of error probability P t ) 
is provided by the difference of the upper-bound and lower-bound 
values of 2E/N P (in dB) for given values of and P. For P. = 10 -4 
and P = 4, the sharpness is \\ dB for = 0.05 and 2\ dB for = 0.1. 
This measure of sharpness appears to be relatively insensitive to the 
value of P. An alternate measure would be the difference in error 
probability for a given value of 2E/N P, and this measure is indeed 
markedly sensitive to P. 

In the region of the curves corresponding to high signal-to-noise 
ratio, there is an improvement in error probability associated with 
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LOWER BOUNDS 
= 0.2- 

= 0.1 - 
= 0.05- 



UPPER BOUNDS 

^P-O.Z 
^„-0=O.I 
0=0.05 




P = 2 



5 10 15 20 25 

2E/N P SIGNAL -TO -NOISE RATIO PER PATH IN DECIBELS 

Fig. 1 — Error-probability bounds for widely-orthogonal signaling, P = 2. 

larger P; the curves become straight lines since P. becomes proportional 
to (2E/N P)~ P . However, this improvement is in part attributable 
to choosing 2E/N P, the average per-path signal-to-noise ratio, as 
the abscissa rather than 2E/N , the total signal-to-noise ratio. To 
obtain plots vs 2E/N , one moves the P = 2 n curves to the right 
by 3n dB; then, the improvement with increased P is less dramatic 
in this region of high signal-to-noise ratio. 

4.3 Computing Spectral Bounds 

It has been observed that bounds on the error probability for the 
case of widely orthogonal signals can be obtained from bounds on the 
spectrum of B (m)<r, m = 1, 2. We now give several easily computed 
formulas for these bounds. 

Recall that B(m) is denned to be / dt b*(t, m)h(t, m), a matrix of 
scalar products or a Gram matrix. In general, this is uninformative, 
since a matrix is a Gram matrix if and only if the matrix is positive 
semidefinite. However, we will shortly use the fact that in our case 
the diagonal entries of B(m) are unity because of the normalization. 
Next, note that B(m)<r is similar to a*B(m)a*, a hermitian matrix 
which has real roots (since a is a real diagonal matrix with positive 
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entries, the matrices a* and a~* exist; then o-'[.B(m)o-]<r~* = <r*B(m)«r*). 
When #(m) is diagonal or nearly so, the roots of B(m)a should be 
close to the entries of a; this is justified by the following theorem. 9 
The characteristic roots of any matrix A lie in the closed region of 
the z-plane consisting of all the disks [z:\z — A it | ^ 2Z,v< | A ti \, 
% =■ 1, 2, ■ • • , P). In our case, the region must be on the real line, 
and we obtain a set of not necessarily nonoverlapping intervals centered 
about {o-,|, the half-widths being {^iv* I B if (m) | a A when we take 
A = B(m)<r. The spectral bounds are then the rightmost right-end point 

max [An + E I A lt |], 

i i f i 

and the leftmost left-end point 

mm[A« - E|4, |] 

(when it is positive). 

A family of spectral bounds is obtainable from this theorem by apply- 
ing it to B(m)a and to matrices similar to B(m)a, e.g., o*B(m)c*, oB(m), 
and more generally <x a B{m)o l ~ a , £ a ^ 1. Thus, we have the family 
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Fig. 2 — Error-probability bounds for widely-orthogonal signaling, P = 4. 
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Fig. 3 — Error-probability bounds for widely-orthogonal signaling, P = 8. 



of upper spectral bounds 

max [<r { + £ff" | B u (m) \ <r)~ 



< a ^ 1. 



(6) 



The question arises: which is the smallest upper bound? It is not true 
in general that a bound is attained for the value of i that maximizes 
o-, , but suppose this is the case when a = 0. That is to say, suppose 
<Ti = max t a k and that 



r.(l+ Z\B< 

L i*i 



(*) 



= max a k 

k 



[l+ E|B w (m)|4 



Then it follows that this is the smallest bound in the family, for ojoi ^ 1 
implies that 



E|£„(m)|^ E|S„("0. , 7 



a ; 



and hence 



r,(l + Z | B w (m) | *] ^ <r,[l + E | fi„(ro) | (fj 

^max^l+gl^MlQ 1 " ]}- 
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Similarly, we have the family of lower spectral bounds 

min {<r, - J2<t- \ £„(m) | a)'"}, Ogagl, (7) 

The largest lower bound is obtained when a = 1 provided that o\ = 
min* cr fc and that 



*.[1 - E \B tl (m) |] = min {<r t [l - £ I *«( 
To see this, observe that o-,/ ". = 1 implies 

Z\B i Am)\fc) 1 ~'' > Z|5,,(m)|, 



///) 



and hence 



*.-[! " E I *■■/(»») |] ^ «rll - D \B,,(m) | (^Y~ j 

;X. L i*i \o-.7 J 

g, r{+-S |B " ( ' K)l fe)"l' 

It should be noted that less sharp bounds are easily obtained. For 
example, the matrix <rB(m) yields the upper bound 

max {«r t -[l + £ I B u {m) |]J ^ max tr,- max [1+21 5 w(™) |], 

and the right-hand side is easily computed. The corresponding lower 
bound is 

min [<r t [l - Z I B „(w) |]j ^ [min <r,][l - max £ | B ti (m) |]. 

These less sharp bounds are easier to compute than those obtained in 
a similar fashion from B(m)a or from a a B(m)<T L ~ a . 

Also, it should be noted that sharper bounds can be obtained by 
employing a sharper theorem of matrix theory: 9 The characteristic 
roots of any matrix .4. lie in the closed region of the z-plane consisting 
of all the ovals | z - A u \ \ z - A„ | ^ (£, k ^ A «)£»-* An), i *■ j- 
We do not pursue these bounds, but note that simple formulas are 
obtained only when all paths have equal strength, o\ = 1/P. 

It is now clear that when B(m) is essentially diagonal, with 
Z^ivt | Bn(m) | « 1 for all i, the path gains o-,- are good nominal values 
for the characteristic roots of B(m)a. If, moreover, these path gains 
are equal, or approximately equal, then the upper and lower spectral 
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bounds are close to one another. When this narrow-spectrum condition 
prevails, the positive and negative portions of the spectrum of L(m)Q 
are also narrow, and the bounds on error probability are sharp. 

V. OTHER SPECIAL CASES 

5.1 Well-resolved Paths and the Theory of Diversity 

We consider the case in which the signals are resolvable, B(l) — 
B(2) = I, i.e., the paths are well separated in time and frequency 
so that any time-Doppler shifted version of a signal is orthogonal to 
any other version of itself. Moreover, we also assume that 5(1, 2) = 
/ dt b*(t, l)h(t, 2) becomes a diagonal matrix, B{1, 2) = pi where 
p = / dt x?(t)x a (t), i.e., the paths are sufficiently separated so that 
any version of one signal is orthogonal to all but the same-path version 
of the other signal. 

It is then easily seen that the covariance matrix is comprised of 
diagonal submatrices. For m = 1, 

L"(l) = c + {NJ2E,)I L l2 (l) = p[a + (JVo/271;,)/] 

L 2, (l) = P *[<r + (N /2E 1 )I] L 22 (l) = | p |* c + WJTBdI- 

For m = 2, assuming E 2 > 0, 

L u (2) = (E 2 /E 1 )[\ P | 2 a + (N /2E 2 )I] 

V\2) = p(E 2 /E,)[<x + {N Q /2E 2 )I] 

L 2, (2) = P *(£ 2 /£i)[" + (N /2E 2 )I] 

L 22 (2) = (E 2 /E,)[a + (N /2E 2 )I]. 

Moreover, the matrix Q is diagonal, being related to 

(2E m /N )H(m) = (2E m /No)[(2EJN )I + a' 1 ]- 1 = a[a + (N /2E m )iy\ 

It then follows that L(m)Q is comprised of diagonal submatrices. 
To find the spectrum, the order of the determinantal equation can be 
reduced from 2P to P. Then the argument of the determinant is quad- 
ratic in X. For the case E, = E 2 , a method of Turin [(22)-(23) in Ref. 5] 
can be used relating the \ k to the eigenvalues (elements) of <r. 

The above example brings the present analysis in contact with the 
theory of diversity combining, see e.g., Ref. 3, Sec. 7.4. Turin, 5 for 
example, considered the case in which separate waveforms are available 
and the fading is nonindependent in general. In our analysis, only one 
waveform is in general available. But in the case of well-separated paths, 
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we may assume P separate signal waveforms have been observed. 
However, these separate waveforms must fade independently in keeping 
with our general discrete-path model, and the on-diagonal component 
matrices of L(m), viz., L ll (m) = {N Q /2E l )(Z l Z\) m and L 22 (m) - 
(N /2E 1 )(Z 2 Zl) m , are themselves diagonal matrices. It is still entirely 
possible that L i2 (?n), the off-diagonal component matrix of L(m), is 
not a diagonal matrix; e.g., when x 2 (t) is a delayed version of Xi(t), 
then time overlap may preclude B(l, 2) being diagonal even though 
J5(l) and B{2) are diagonal. But when we assume that 5(1, 2) is also 
diagonal, then we obtain the form for L(m) exhibited above. It can be 
observed that this is precisely the result Turin obtained for the case 
of optimum diversity combining, where his not necessarily diagonal A 
becomes our diagaonal a. When B(l, 2) is not diagonal, then our 
results do not specialize to the form given by Turin, a reflection of 
the fact that the multipath channel is not in general fully equivalent 
to a diversity channel. 

5.2 On-Off Keying 

Another example is the case of on-off keying in which E 2 = 0. The 
test statistic Z*QZ becomes 

[(N /2E i ) i Z(l)]\2E l /N )H(l)[(N /2E 1 ) i Z(l)], since Q 22 = 0. 

Thus, the distribution is determined by the spectrum of the matrix 
L"(m)Q u , where 

L n O) = 6 ml B(l, m)aB(m, 1) + (N /2E 1 )B(l) l 
Q n = ZT'(l)!/ + WofiBdUHWT 1 . 

Observe that we no longer have the difference of positive-definite 
forms, the test statistic now being a positive random variable. The 
threshold (N /2E l )d(2) is 



(ATo/2/s.) log det 



(f)B(l). + / 



.ViVo 



which is positive since the eigenvalues of {2EJN Q )B{\)a -f / are 
greater than unity. 

Assuming that the spectrum of V\m){2E l /N )H{\) lies in the 
interval (m, m), where n and y. are functions of m, the bounds on the dis- 
tribution function are 



am-\N /2E)e;P,0] ^ F n 



■IE 



^ G[b)-\N /2E)d;P,0]. 
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Recall that G(x; P, 0) is related to the incomplete gamma function, 
G(x;P,0) = I(x,P - 1). 

The spectral bounds must exhibit two forms of (2fi 1 /iV )-dependence. 
When m - 1, L"(l)Q n = B(l)a, and bounds on B(l)a become m and ji. 
When m = 2, L"(2)Q" = (N /2E X )[I + (JWaEOPttM" 1 }" 1 . so that 



\-ii-i 



-\-ii-i 



p. = (JSTo/2^0[l + {No/2E l )(o,)' 1 ] 

where the spectrum of B(l)a is confined to (co, d>). 
Collecting our results, when m = 1, 

F 1 [(N /2E)d] ^ I{( iS )-\No/2E)P log [(2E/N Q )Z> + 1];P - 1) 

F 1 [(N /2E)d] ^ H(u)-\No/2E)P log [(2E/N )c, + 1];P - 1}. 

Similarly, when wi = 2 

F 2 [(N /2E)d] ^ 7{[1 + (N /2E)(aT 1 ]P log [(2E/N )a + 1];P - 1} 

F,[{No/2E)0\ ^ H\X + (No/2E)(o>T 1 ]P log [(2#/tf )g + 1];P - 1). 

These results permit the computation of error-probability-bound curves 
that would be universal in the same sense as the curves for widely- 
orthogonal signaling, i.e., the curves would apply to any element of 
the set of channels for which the spectral bounds are met. 

VI. CHERNOFF BOUNDS 

6.1 General Case 

Up to this point, consideration of spectral bounds has lead to error- 
probability bounds which are sharp when the spectrum comprises 
narrow positive and negative portions. These bounds are easy to employ 
when 5(1, 2) = and B(l), B(2) are nearly diagonal matrices. But 
in more general cases, the estimation of spectral bounds may be difficult 
and bounds may be poor approximations of eigenvalues. We turn to 
another technique of bounding error probability which does not ex- 
plicitly require spectral bounds. 

Consider the error probability when hypothesis m = 2 is true, 
P,(2) = Pr {Z^QZ > (N /2E)d | 2}. Recall that the unit step function 
U(x) is unity for x > 0, zero for x < 0, and one-half for x = 0. Then 
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P.(2) = Pr [U[Z*QZ - (N /2E)6] = 1 | 2} 

= Z 2 \V\Z'QZ- (N /2E)d]}, 

where 8 2 denotes expectation under hypothesis m = 2. But since 
U(x) ^ exp (/i 2 .r) for any /z 2 > 0, we have 

P.(2) ^ £ 2 {exp M2 [^ t (?2 - (^ /2£)fl]}. 

This average can readily be computed, since Z^QZ has the same dis- 
tribution as X A *( 2 ) I z t \ 2 > where (Xi(2)} is the spectrum of L(2)Q. 
Since 8,z k = 0, 82,2? = 8, k , and Sz/3* = 0, the Gaussian variables 
{Re Zj\, {Im 2,} are independent with zero mean and variance equal 
to 5. Thus, P e (2) is bounded from above by 

exp [- M8 (#o/2£)0] IT £ cxp ( M2 X,(2) | Re 2* | 2 ) , 

_A = 1 J 

where the outer square appears because the product involving {Im z k \ 
has been suppressed. But a standard calculation shows, 

8 exp ( M2 A,(2) I Re z k | 2 ) = [1 - MiX*(2)]"*i when M2 X,(2) < 1, 
and our bound is 

exp [-p 2 (N o /2E)0] il [1 - M 2 X t -(2)]-'. 

Thus, 

P.(2) ^ exp [- M2 (^ O /2E)0] dot" 1 [I - m 2 L(2)Q], (8) 

which holds for all ps such that < n» < [max A,, (2) ]-*. 

The above procedure is adopted from the technique due to Cher- 
noff (see Rcf. 3, Sec. 2.5 and 7.4). Here, we do not have identically 
distributed variables; indeed, half are positive and half are negative 
random variables. 

To find the best value of ix->, we write the bound as 

expj-^0- In II [1 — ixM?i\ 

and differentiate the argument of the exponential. A necessary con- 
dition for an extremum is that the derivative be zero, and this yields 



= E 7-t- = to {[{L{2)Qr - pJY 

k „l Aftdl — H2 



1622 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1967 

If the value of /x 2 that satisfies this equation lies within the allow- 
able interval [0, mur*A*(2)] f then this value of /* 2 minimizes the 
upper bound. A minimum occurs because the second derivative of the 
argument of the exponential is positive, being 

f r x*(2) i 3 

fe U - M 2 X*(2)_ " 
In a similar fashion, the error probability for m = 1 can be over- 
bounded. 

P.(l) = Pr {-Z'QZ > -(N /2E)d | 1} 

fg S,{exp n x [-Z'QZ + (N /2E)6]\ 

P.(l) ^ exp MN O /2E)0\ det" 1 [7 + ^,L(1)Q]. 
The best value of /xi satisfies 

(N /2E)9 = tr {LQ)Q[I + nMVQT 1 }, 
provided this value lies in the allowable interval [0, max _1 (— A fc (l))]. 

6.2 Widely-orthogonal Signals 

Consider the case in which the signals are widely orthogonal, 
P(l, 2) = 0, but have equal energy, E y = E 2 = E, and arc equilikely, 
a} = a2 = i. The overbound on P r (1) is obtained from the spectrum 
of L(1)Q which comprises the spectrum of L 11 (DQ 11 together with 
the spectrum of L 22 (1)Q 22 . Thus, 

P.O.) ^ exp MN /2E)ff] de^ 1 [I + ^"(l^def [7 + ^L 22 (l)Q 22 ]. 
But the matrices used here were related in Paragraph 4.1 to P(l)o- and 
B(2)a, and our bound becomes 

exp [ Ml (AT„/2E)0] det" 1 [7 + ^B(l>] 

•det- {7 - ^(N /2E)[I + {Nj2E)(B(2)<rVr\. 
After some manipulation, this bound becomes 

det [p(2)<t + (§) 7] 
'*{*(£)*». + ($ )/]«{£ - „(|)>(2), + (f ) 
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where 

exo \u (N /2E)9] - <f d0t \B 8* + W t/mA" "™ 

exp MN n /2L)e\ - ^ et [B{2)a + {No/2E)I] j 

The maximum allowable value of p.\ is determined by the largest 
eigenvalue of L 2Z (1)Q 22 which in turn is determined by the largest 
eigenvalue of B(2)a: 

< m, < y + max_1 (**)> 

where {8 A } is the spectrum of 7i(2). 

The best value of m is found from the relation 

(N Q /2E)6 = 2 - , X * (1 j> 

-tr{L"(DQ ,1 [J + MiL ,, (l)0 Ii r) 

+ tr iL 22 (l)Q"[/ + Ml L 22 (l)Q ,, ]- , i, 

where we again have exploited the decomposition of the spectrum of 
L{1)Q. After some manipulation, we find 

(N /2E)d = tr \B(l)a[I + M.fiOVr 1 ) 

- tr {fi(2)a[j + (^ - Mi )b(2).]"'}. 

An approximate solution can be obtained for the case of high signal- 
to-noise ratio. Let /zj = /Z,(2£ , /A r y ); the relation becomes 

(N /2E)d = Z i + {2E/N )fiv k ™ ^ 1 + (2B/Jvi(l - *)«»" 

Suppose jB,(2#/JV )a> 4 » 1 and (1 - fi l )(2E/N )8 k » 1. Then the 
right side becomes approximately 



(2E/N ) Pl (2E/N )(1 - ft) 

Equating this to (N /2E)6 and solving the resulting quadratic for the 
root applicable for the case = yields 



Mi = 



1 + &+V+4P 
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When 0/2P is small, this value of ji x is approximately 

2 L 4PJ ' 

and the corresponding value of Mi is (E/N )[l - (d/P)] which is ap- 
proximately at the midpoint of the allowable interval. 
In a similar fashion, the overbound on P,(2) is 

exp [-y. 2 {N /2E)d] det" 1 [I - n 2 L 11 (2)Q u ] det" 1 [I - /i 2 L M (2)0 M L 

which becomes 

exp [-» 2 (,N /2E)6] det"' [I - n 2 (N /2E)[I + (N /2E)(B(l)af l T 1 } 

•det" 1 [7 + n 2 B(2)<r], 
or 

E \- p { ( No 



det [*(!), + gj) /] 



■--Kl--.fi)B(«> + &)/}*[^a)^ + *)I] , 

where 

_ J det[^(2) t 7 + (iVn/2E)J] V- (Ar ° /2g) 



exp [- M ,(AT o /2£)0] = \ det [5(1)ff + (tf o/ 2W 

The maximum allowable value of y. 2 is determined by the largest eigen- 
value of L ll (2)Q n in turn determined by the largest eigenvalue of 
B(l)a: 

< M2 < tt + max -1 (w*), 

where {o) fc } is the spectrum of B(1)<t. The best value of n 2 satisfies 
(N /2E)6 = tr fL n (2)Q n [7 - M2 L 1, (2)Q n ]- 1 } 

+ tr{L 22 (2)Q 22 [7- M2 L"(2)Q 2 T 1 } 

(AT o /2£)0 = tr |fl(l)«r (^ - M»)tf(l)<r + 7 J } 

-tr {B(2)a[I + t x 2 B(2)aT 1 }. 
Let a 2 = p 2 (2E/ N ) and suppose p 2 (2E/N )8 k » 1, (l-/B 8 )(2#/tf K» 1. 
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Then the right side becomes 

P 



(2E/N )(1 - fc) (2E/N )n 2 ' 
and the approximation of the best value of /Z 2 is 



P 2 = 



i-^J-U- 



4P 2 



which is approximately |[1 - (0/4P)] when (0/2P) « 1. 

The foregoing results can be specialized to the case in which the 
paths are resolvable, 5(1) = B(2) = I. Then = 0, and it is easily 
seen that the best value of /Z m is h. Both overbounds become 

(F/ o N s-e dot [a + (N /2E)I] 

and this agrees with equation 7.134 in Ref. 3. 

It should be noted that /Z m = h is always an allowed value of p m . 
For the case of resolvable paths, it is the best value, and whenever 
0/P « 1 and 2E/N is sufficiently large, it is close to the best value. 
Using p m = h, we can obtain an overbound for both error probabilities, 
i.e., for P.(m), m = 1, 2. This overbound is 

(E/2NX" exn ft I \) det [B(3 - m)v + (N /2E)I] 

' ' [ K2 ' V det [B(l)a + {N Q /E)I} det [B{2)a 4- {N Q /E)T\ 

(9a) 

The factor exp Q | |) can also be written in terms of determinants. 
When det [£(l)<r + (N /2E)I\ is larger than det [B(2)<r + (N /2E)I], 
we have 

fdet [£(1V 4- (N a /2E)r\ \* 

Idet [ff(2)o- + (N /2E)r\f ' 

and when the reverse inequality holds, exp (J | |) is the reciprocal 
of the above. 

For the case in which the spectrum of B(m)a lies in the interval 
(1-/3/P, 1+jl/P) where p < 1, the overbound can be further over- 
bounded. The factor involving determinants is less than 



exp (i | |) = 



(9b) 



t+'+ftF) 



b - ' + <¥) 
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and | 8 |/2 is less than 

P 1 + + (N P/2E) 
2 10g 1 - + (N P/2E) 

It follows that the Chernoff bound is less than 

4[1 + + (NoP/2E)] 



(N P/2E)' 



(10) 



,[1 - + (iV P/2^)Hl - IS + 2(AT P/2^)] 2 
Numerical values of this bound are given in Fig. 4, and it has the 
same general character as the spectral-related bounds. Rather than 
sharpness given a nominal value of error probability P e , we consider 
the sensitivity measured by the change in 2E/N P (in dB) vs /?; 
for P e = 10" 4 and P = 4, the sensitivity is 2 dB for p = 0.1. The 
sensitivity does not markedly increase with an increase in P, in agree- 
ment with the behavior of the sharpness of the previous bounds. 

Comparison of the Chernoff bound with the previous bounds is 
conveniently done for the case = (cf. Sec. 7.4 of Ref. 3). The 
Chernoff bound does not specify a signal-to-noise ratio (required to 
achieve a nominal P e ) excessively greater than the previous value; 
for P = 4, less than 2.2 dB difference is observed. This excess does 
decrease with increasing P. Moreover, it is entirely conceivable that in 
a broad-spectrum case with a large number of paths, an exact value of 
the Chernoff bound would be better than the spectral-bound result. 
Of course, our inexact (overbounded) Chernoff bound is poor in the 
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Fig. 4 — Overbounded Chernoff bounds for widely-orthogonal signaling. 
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broad-spectrum case, but a Chernoff bound using the proper values 
of the determinants should be good for two reasons, (i) Such a bound 
reflects the precise values of the eigenvalues of the matrix L(m)Q. 
(it) When P is large, the probability density function is bell shaped 
with the probability ''mass" being concentrated near the mean and 
most of the tail mass being at the leading portion of the tail; then 
the tail mass can be weighted by the exponential function with little 
error. On the other hand, the spectral-bound approach suffers in the 
broad-band case since the spectral bounds are not meaningful approxi- 
mations of all the eigenvalues. 



VII. DISCUSSION 



Having observed that exact computation of error probability is 
cumbersome and depends upon an often inordinately large number of 
parameters, we considered error-probability bounds (2) that are uni- 
versal in the sense that they apply to any one of a set of channels 
satisfying spectral bounds (1). Our bounds employ (3), the distribu- 
tion function of the difference of chi-square variables. For the special 
case of widely orthogonal signals, we obtained bounds employing 
parameters (5) in terms of the spectral width p, see (4), of the 
matrices B(m)<r. Plots of these bounds showed that sharpness meas- 
ured in dB change of 2E/X„P with respect to (3 for a fixed value of 
error probability is not sensitive to the value of P. We presented a 
technique for obtaining spectral bounds for B(m)a when it is nearly 
diagonal, representative results being (6.) and (7). This technique can 
also be applied to L(m)Q for the more general case in which the 
signals are not widely orthogonal. 

The case of resolvable signals {B(m) = I) made contact with the 
theory of diversity; we found that for the multipath channel to be a 
diversity channel, £(1,2) must also be a diagonal matrix. Of course, 
the previous results also were in contact with diversity theory. With 
B{1, 2) = (a diagonal matrix) but B{m) not necessarily diagonal, 
our results generalize those of diversity theory in the following sense. 
The special case fi = corresponds to a diversity channel with equal 
link gains, but the general case p ¥= can arise in the nondiversity 
situation when the matrix B(m) is not diagonal. (If B{m) were 
diagonal, B(m) = I and the diversity case prevails.) 

We then turned to the Chernoff bound (8) which does not ex- 
plicitly employ spectral bounds. The overbounded form (10) for the 
case of widely orthogonal signals was poorer than the previous bound 
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when (3 = 0. Nevertheless, there is promise that in a broad-spectrum 
case, the original form (9) would be better than the spectral-related 
bounds. A further advantage is that once the determinants are eval- 
uated, perhaps on an electronic computer, the error-probability bound 
is immediately obtained. In contrast, the spectral-related bounds 
require a certain amount of computation involving incomplete gamma 
functions even after spectral bounds are obtained. 
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APPENDIX A 

Here we show that the number of positive eigenvalues of LQ 
equals the number of negative eigenvalues.* Recall that L is positive 
definite and that Q can be written in the partitioned form 

fa" 



q = 



where Q n and — Q 2Z are positive definite. Clearly, the number of positive 
eigenvalues of Q equals the number of negative eigenvalues. We can 
construct a family of positive definite matrices L t , < t £ 1, such 
that L = I, Li = L, and L t is continuous in t. For example, let L, = 
(1 — t)I -f- tL; L, has positive eigenvalues {(1 — t) + ty k \, where 
[y k ] are the eigenvalues of L. Now the eigenvalues of L,Q are real, 
for L,Q is similar to the Hermitian matrix L\QL\ = Lj^L^L^, where 
L\ and L~^ exist since L t is positive definite, Moreover, the eigenvalues 
of L t Q are continuous in t, since L, is continuous in t. But L,Q never 
has a zero eigenvalue, for L, is positive definite and (L t Q)~ l = Q~ l Lj* 
always exists. Since the eigenvalues are real, continuous in t, and never 
zero, it follows that no positive eigenvalue of L Q can become negative 
as t varies on [0, 1], and no negative eigenvalue of L Q can become 
positive. The conclusion is established. 

APPENDIX B 

This appendix presents another derivation of the distribution func- 
tion of 2^f J z k | z — a ]Cp+! ! z k | 2 - This derivation makes contact with 
* We are indebted to B. H. Bharucha for the conception of this proo£ 
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the special functions that have appeared in analyses of diversity chan- 
nels; also, this derivation appears to admit generalization to the case 
z k = Re z k with (ZjZ k ) = 8 ik . (An odd number of variables in the real 
case corresponds to half-integer P in the complex case.) 
The density function of 2jf I ** | 2 ' s 



{x > 0), 

(x < 0), 

h I is 

(x > 0), 

(x < 0). 



K 


-v 


r 


-1)! 


. 






and the density function of — a /.V* , 





g(x) = 



(-•r) J 



la'(P- 1)! 
The density of the sum is the convolution of the densities, 



h(x) = I dy j(,y)g(x - y), 



where the first argument of max ( • , • ) arises from the truncated form 
of / and the second argument arises from the truncated form of g. 
It follows that 

= exp(.r/a) r _ r-> r_ (1 + A 1 

CX r (P - 1)!(P - l)!4-x<0.x) " _ V« / _ 

For the case x > 0, the lower limit is x. For the case x < 0, the integral 
can be cast into the form of the integral for the case x > by a change 
of variable. The result differs only in the exponential factor, i.e., 



h(x) = 



exp (— x) 



a p (P - 1)!(P - 1)! 



[jy{y- |.r|) / '-V'- | exp[-Q+l) ; /], 



x < 0. 



The integral can be evaluated with the aid of relation (12) on page 202, 
Vol. II of Ref. 10, and the common result for the cases x < and x > is 



x \ P * exp 



h(x) = —j= 



Ml 



~ Kp-i 



1 + a I x 



VVa(l +a)''-KP - 1)! 
where K P -\(z) is the modified Bessel function of the third kind. 
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The above expression for the density is valid for all P, noninteger 
as well as integer. But in our application, P is an integer; a relation 
on page SO of Ref. 11 yields 



The density is then 
h{x) = fe) GXP 



2z fe*!(P- l-fc)!(2*)* 
1 — a \ X 1 + a I x 



•> 



" fa (P - l)!fc!(P - 1 - k)\ \1+ «/ a 

When x < 0, the exponential becomes exp (s/a), and when x > 0, 
it becomes exp ( — .t). 

Observe that when a = 1, the density is symmetric. When a < 1, 
the factor exp [(1 — a/a)x/2] shifts the mass to the right. When a — » 0, 
it can be shown that /i(a;) — * /(x). 

To obtain the distribution function G(y; P, a), consider first the 
case y < 0. Since /!«, dx h(x) equals /"„, da: h(-x), the following integral 
arises in each term of the sum, 

f fe.-^-feV" M - (P - 1 - fc)![l - 7(| y \/a, P-\- k)]. 

The case y > is treated by considering J° M rfx- fe(a;) + Jl dx h(x). 
The integral that arises is just (P - 1 - fc) U(y, P - I - k). These 
steps establish our final result, quoted above. 

Our result could also have been obtained from the Fourier transform 
of the characteristic function (1 - fl)" p (l + tte)"'. The Fourier 
transform of (a + it)~ 2 "(fi - it)' 2 " is given by relation (12) on page 119, 
Vol. / of Ref. 10 in terms of Whittaker functions that reduce to Bessel 
functions for the case ju = v = P/2 in view of relation (14) on page 265, 
Ref. 12. The density function can thus be obtained. 
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